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i i ■ Supersymmetric models with bilinear ii-parity violation (BRPV) can account for the observed 

neutrino masses and mixing parameters indicated by neutrino oscillation data. We consider minimal 
supergravity versions of BRPV where the lightest supersymmetric particle (LSP) is a neutralino. 
, This is unstable, with a large enough decay length to be detected at the CERN Large Hadron 

Collider (LHC). We analyze the LHC potential to determine the LSP properties, such as mass, 
lifetime and branching ratios, and discuss their relation to neutrino properties. 
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• rH ' I. INTRODUCTION 

X: 

Elucidating the electroweak breaking sector of the Standard Model (SM) constitutes a major challenge for the 
Large Hadron Collider (LHC) at CERN. Supcrsymmetry provides an elegant way to stabilize the Higgs boson scalar 
mass against quantum corrections provided supersymmetric states are not too heavy, with some of them expected 
within reach for the LHC. Searches for supersymmetric particles constitute a major item in the LHC agenda [ll4l0|. 
as many expect signs of supcrsymmetry (SUSY) to be just around the corner. However the first searches up to ~ 
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5 fb 1 at the LHC interpreted within specific frameworks, such as Constrained Minimal Supersymmetric Standard 
Model (CMSSM) or minimal supergravity (mSUGRA) indicate that squark and gluino masses are in excess of ~ 1 



TeV 



Hi. 



Despite intense efforts over more than thirty years, little is known from first principles about how exactly to realize 
or break supersymmctry. As a result one should keep an open mind as to which theoretical framework is realized in 
nature, if any. Supersymmetry search strategies must be correspondingly re-designed if, for example, supersymmetry 
is realized in the absence of a conserved R parity 0, [T3| • 

Another major drawback of the Standard Model is its failure to account for neutrino oscillations 0, , whose 
discovery constitutes one of the major advances in particle physics of the last decade. An important observation is 
that, if supersymmetry is re alized w ithout a conserved R parity, the origin of neutrino masses and mixing may be 
intrinsically supersymmetric 
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Indeed an attractive dynamical way to generate neutrino mass at the weak scale is through non-zero vacuum 

THE I 



This leads to the minimal effective 



expectation values of SU(3) <B> SU(2) Cg> U(l) singlet scalar neutrinos 
description of R parity violation, namely BRPV 22 1 . In contrast to the simplest variants of the seesaw mechanism [2 

Here we 



24 



such supersymmetric alternative has the merit of being testable in collider experiments, like the LHC 
analyze the LHC potential to determine the lightest neutralino properties such as mass, decay length and branching 
ratios, and discuss their relation to neutrino properties. 



II. BILINEAR R-PARITY VIOLATING SUSY MODELS 

The bilinear R-Parity violating models are characterized by two properties: first the usual MSSM R-conserving 
superpotential is enlarged according to [E| 

WbRPV = WmSSM + SabtiLiHt > (1) 

where there are 3 new superpotential parameters (ej), one for each fermion generation . The second modification is 
the addition of an extra soft term 

Ksoft = Vmssm — SabBieiL^H^ (2) 

that depends on three soft mass parameters B^. For the sake of simplicity we considered the R-conserving soft terms 
as in minimal supergravity (mSUGRA). Notice that the presence of the new soft interactions prevents the new bilinear 



2S 



terms in Eq. ([T]) to be rotated away 

The new bilinear terms break explicitly R parity as well as lepton number and induce non-zero vacuum expectation 
values Vi for the sncutrinos. As a result, neutrinos and ncutralinos mix at tree level giving rise to one tree-level 
neutrino mass scale, which we identify with the atmospheric scale. The other two neutrino masses are generated 



through loop diagrams 3l|, |32J . This model provides a good description of the observed neutrino oscillation data 14 



The BRPV-mSUGRA model is defined by eleven parameters 

m , TO1/2 , tan/3 , sign(^) , A , e; , and B t , (3) 

where m 1 / 2 and mo are the common gaugino mass and scalar soft SUSY breaking masses at the unification scale, Aq 
is the common trilinear term, and tan/3 is the ratio between the Higgs field vacuum expectation values (vevs). In our 
analyzes the new parameters (ej and Bi) are determined by the neutrino masses and mixings, therefore, we have only 



1 In a way similar to the /i term in the MSSM superpotential, the required smallness of the bilinear parameters ej could arise dynamically, 
through a nonzero vev, as in [l9H2ll . [29ll and/or be generated radiatively [30(1 . 
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Figure 1: Lightest neutralino decay length as a function of mSUGRA parameters mo and m 1 ^ 2 , for Ao 
and sign(/x) > 0. 



-100 GeV, tan/3 = 10 



to vary the usual mSUGRA parameters. For the sake of simplicity in what follows we fix Aq = —100 GeV, tan/3 = 10 
and sign(^) > and present our results in the plane mo ®m[/ 2 . 

Due to the smallness of the neutrino masses, the BRPV interactions turn out to be rather feeble, consequently the 
LSP has a lifetime long enough that its decay appears as a displaced vertex. We show in figure [1] the LSP decay 
length as a function of too and mi/2, when the remaining values for sign(//), A and tan/3 are taken as mentioned 
above. Therefore, we can anticipate that the LSP decay vertex can be observed at the LHC within a large fraction of 
the parameter space. 

Depending on the SUSY spectrum the lightest neutralino decay channels include fully leptonic decays 



X? 



and 



x? 



with £ = e or fi; as well as semi-leptonic decay modes 



X? 



vqq 



Xi -> Tq'q 



X?^ 



' q and Xi ~~ vbb 



If kinematically allowed, some of these modes take place via two-body decays, like y£ 



~° 1 Zv, or x°i 



Yi 
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W^y^, Xi ~> W^t*, 
In addition to these 



hu, followed by the Z, or h decay; for further details see Ref. 
channels there is also the possibility of the neutralino decaying invisibly into three neutrinos, however, this channel 
reaches at most a few per-cent [j^ 2 . 

Neutrino masses and mixings as well as LSP decay properties are determined by the same interactions, therefore, 
there are connections between high energy LSP physics at the LHC and neutrino oscillation physics. For instance, 
the ratio between charged current decays 



Br(x; -» W ± ^) 



(4) 



of 6*23 and Ato| 2 from Ref. 39]. 



is directly related to the atmospheric mixing angle [38( , as illustrated in the right panel of figure [2] this relation was 
already considered in Ref. 27 L The vertical bands in figure[5]correspond to the latest 2a precision in the determination 



2 However, in models where a Majoron is present, it can be dominant 134437(1 
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Figure 2: Correlating LSP decay properties with neutrino oscillation parameters. The left panel shows the connection between 
the displayed LSP decay length parameter and the atmospheric squared mass scale Am§2- The right panel depicts the relation 
between Br(x? — ► W ± ^t T )/Br(xi — ► W /± r =F ) and the atmospheric mixing angle. The vertical shaded bands indicate the 2a 
allowed values of the corresponding neutrino oscillation parameters [3^ . 



Another interesting interconnection between LSP properties and neutrino properties is the direct relation between 
neutrino mass squared difference Am§ 2 and the ratio 

Br(x° -> W£ + Zv) 

as is illustrated in left panel of figure [2] Here Lq is the LSP decay length and one has to sum over all leptons and 
neutrinos in the final states. One can understand this relation in the following way. In the BRPV model the tree-level 
neutrino mass is proportional to rn^ rec oc |A| 2 , where |A| 2 = ^ A 2 , with Aj = evd + (iVi, is the so-called alignment 



^32 - 57730 ; — e ; (5) 



vector. Couplings between the gauginos and gauge bosons plus leptons/neutrinos are proportional to Aj as well 33 1. 
Thus, one expects that after summing over the lepton generations the partial width of the neutralino into gauge 
bosons is also proportional to |A| 2 . The decay length is the inverse width and dividing by the branching ratio into 
gauge boson final states picks out the partial width of the neutralino into gauge bosons. This leads to the correlation 
of i?32 with the atmospheric neutrino mass scale, since m-Atm is identified mostly with m^ 100 , apart from some minor 
1-loop corrections. 



III. ANALYSES FRAMEWORK AND BASIC CUTS 



Our analyzes aim to study the LHC potential to probe the LSP properties exp loring its detached vertex signature. 

where all the properties of our 



exp 
4fl 
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Wc simulated the SUSY particle production using PYTHIA version 6.408 
BRPV-mSUGRA model were included using the SLHA format 1421. The relevant masses, mixings, branching ratios, 
and decay lengths were generated using the SPHENO code 43J, |44j . 

In our studies we used a toy calorimeter roughly inspired by the actual LHC detectors. We assumed that the 
calorimeter coverage is \r)\ < 5 and that its segmentation is Arj <E> Aip = 0.10 x 0.098. The calorimeter resolution was 
included by smearing the jet energies with an error 

AE 0.50 
~E~ ~~ 



E 



®0.03 



Jets were reconstructed using the cone algorithm in the subroutine PYCELL with AR 

tLn = 2 GeV. 



0.4 and jet seed with a 



minimum transverse energy E^ l n 
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Our analyzes start by selecting events that pass some typical triggers employed by the ATLAS/CMS collaborations, 
i.e. an event to be accepted should fulfill at least one of the following requirements: 

• the event contains one electron or photon with px > 20 GcV; 

• the event has an isolated muon with pr > 6 GeV; 

• the event exhibits two isolated electrons or photons with pt > 15 GeV; 

• the event has one jet with transverse momentum in excess of 100 GeV; 

• the events possesses missing transverse energy greater than 100 GeV. 

We then require the existence of, at least, one displaced vertex that is more than 5tr away from the primary 
vertex 25j - that is, the detached vertex is outside the ellipsoid 

5S X y) + (s^y) + (54) ' ^ 
where the z-axis is along the beam direction. We used the ATLAS expected resolutions in the transverse plane 
(S xy = 20 /im) and in the beam direction (S z = 500 /im). To ensure a good reconstruction of the displaced vertex we 
further required that the LSP decays within the tracking system Lc^ within a radius of 550 mm and z-axis length of 
3000 mm. In our model the decay lengths are such that this last requirement is almost automatically satisfied; see 
figure [TJ 



IV. LSP MASS MEASUREMENT 

In order to accurately measure the LSP mass from its decay products we focused our attention on events where 
the LSP decays into a charged lcpton (e^ or /i ± ) and a W that subsequently decays into a pair of jets. In addition 
to the basic cuts described above we further required charged leptons to have 

p e T > 20 GcV and |r#| < 2.5 . (7) 

We demanded the charged lepton to be isolated, La the sum of the transverse energy of the particles in a cone 
AR = 0.3 around the lcpton direction should satisfy 

e t < 5 GcV . (8) 

Ai?<0.3 

We identified the hadronically decaying W requiring that its decay jets are central 

y T >20GeV , \rij\ < 2.5, (9) 
and that their invariant mass is compatible with the W mass: 

1 17,- 1 < 2.5 and \M ]3 - M w \ < 20 GeV . (10) 

In order to obtain the LSP mass, we considered points in the mo <8> m^ plane with more than 10 expected events 
for an integrated luminosity of 100 fb _1 . We have performed a Gaussian fit to the lepton-jet-jet invariant mass; as 
an illustration of the lepton-jet-jet invariant mass spectrum sec figure [3] As we can see from this figure, the actual 
LSP mass (101 GeV) is with 1% of its fitted value (100.4 GeV). 

In order to better appreciate the precision with which the LSP mass can be determined for other choices of mSUGRA 
parameters we have repeated the analysis for a wide grid of values in the mo (£> m^ 2 plane. The left panel of figure 0] 
depicts the achievable precision in the LSP mass measurement for an integrated luminosity of 100 fb -1 as a function 
of mo <£> m.x/2 for Aq = —100 GcV, tan/3 = 10, and sgn(//) > 0. As one can see the LSP mass can be measured with 
an error between 10 and 15 GeV within a sizeable fraction of the (mo ® JTI1/2) plane. Only at high m 1 / 2 there is a 
degradation of the precision due to poor statistics. The right panel in figure [4] shows that indeed this is enough to 
determine the LSP mass to within 5 to 10% in a relatively wide chunk of parameter space. 



6 




100 120 140 160 180 200 



m z (GeV) 



Figure 3: Illustration of the lepton-jet-jet invariant mass spectrum fitted to obtain the LSP mass. In this figure we considered 
m = 250 GeV, m 1/2 = 250 GeV, tan/3 = 10, A = -100 GeV, and sgn(^) > which leads to a LSP mass of 101 GeV. 

V. LSP DECAY LENGTH MEASUREMENT 



Another important feature of the LSP in our BRPV-mSUGRA model is its decay length (lifetime). Within the 
simplest mS UGRA bilinear R parity violating scheme this is directly related to the squared mass splitting Am| 2 , well 
measured in neutrino oscillation experiments 39( . In this analysis we considered events where the LSP decay contains 
at least three charged tracks, La the LSP decays into £jj, with I = e or /j,. Here we sum over all jets as well as over 
X? — > tW — > Ijj and all three body decays leading to the same final state. 

In figure \5\ we depict the average distance traveled by the LSP as observed in the laboratory frame. As we can 
see, a substantial fraction of the LSP decays takes place within the pixel detector, except for very low toi/2 values. 
It is interesting to notice that the pattern shown in the figure is similar to the one in figure [I] as we could easily 
expect. Since most of the LSP decays occur inside the beam pipe we can anticipate a small backgrounds associated 
to particles scattering in the detector material. 

In order to obtain the LSP decay length (Lq) from the distance traveled in the laboratory frame (d) we considered 
the m b s d/p i) S distribution, with m f, s (p bs) being the measured invariant mass (momentum) associated to the 
displaced vertex, and then we fitted it with an exponential 



where the fitting parameter (Lq) is the LSP decay length. 

In order to disentangle the energy and momentum uncertainties and the statistical errors from the intrinsic limitation 
associated to the tracking we first neglect the latter one. In the left panel of figure [6] we present the expected precision 
in the decay length determination in the plane mo <8> W1/2 f° r an assumed integrated luminosity of 100 fb _1 . As one 
can see, these sources of error have a small impact in the determination of the decay length, except for heavier LSP 
masses where we run out of statistics. In fact, the contribution of these sources of uncertainty is smaller than 5% for 
neutralino masses up to 280 GeV (mi/2 — 700 GeV). 

Clearly, the actual achievable precision of LSP lifetime determination at the LHC experiments depends on the 
ability to measure the LSP traveled distance in the laboratory. We present in the right panel of figure |B] the attainable 
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Figure 4: The left panel presents the error (ctm) on the LSP mass as a function of the mo ® mi/2 point for Aq — —100 GeV, 
tan/3 = 10, sgn/i > and an integrated luminosity of 100 fb~ , while the right panel displays the relative error in the LSP 
mass determination <jm/M^o. 



precision on the decay length assuming a 10% tracking error |45| in the LSP flight distance to get a rough idea. 
Clearly the precision in the decay length gets deteriorated, however, it is still better than 15% within a relatively 
large fraction of the parameter space under this assumption but would get corresponding worse if this uncertainty 
were larger. 



VI. LSP BRANCHING RATIO MEASUREMENTS 

As we have already mentioned, the neutrino mass squared difference Am| 2 controls the ratio given in Eq. [SJ 
therefore we should also study how well the neutralino LSP decay ratio into iW and vZ can be determined. In order 
to illustrate the LHC capabilities in probing LSP properties at high energies we present the reconstruction efficiency 
for the benchmark scenario 



'1/2 



250 GeV and m = 250 GeV 



that yields a rather light LSP (tolsp — 101 GeV) and heavy scalars. For this point in parameter space the LSP 
possesses a decay length cr = 30 /jm and its dominant decay modes have the following branching ratios: 

BR(x? -> W ± e^) = 0.2%, BR(x? -> W ± ^) = 27.6%, BR(x? -> W^) = 31.3%, 
BR(x?-s>&H = 7.1%, B'Rtxl^Zv) = 11.9%, 

BRfx^e^^) = 5.5%, BR(£? -> n^v) = 5.5%, BR(x? -> ^t^v) = 9.5% 

We present in Table U the reconstruction efficiencies of the LSP decay modes for our chosen benchmark point. The 
reconstruction efficiencies for final states containing t's are much smaller, as expected, leading to a loss of statistics 
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in these final states. For an exhaustive study of the reconstruction efficiencies see Ref. 

Wc present in figurc[7]the expected error on the LSP branching ratio Br(x? — > IW + vZ) as a function mo® m i/2 f° r 
an integrated luminosity of 100 fb _1 . In order to evaluate this error we studied the reconstruction efficiency for this 



8 




A A A A • 

AAAAAAAAAAACA* 
AAAAAAAAAAAA«AA 
AAAAAAAAASA****** 
AAAAA*«********** 



■ TTTTTTTTTTTTTTTT 
TTTTTTTTTTTTTTTTT 

\ 

< 0.10 mm ■L lab <1.0mm 

< 0.50 mm ▼ L lab < 5.0 mm 

I i i i I i i i I i i i I i i i I 

200 400 600 800 1000 

m (GeV) 



Figure 5: Average distance traveled by the LSP in the laboratory frame as a function of the mo <8> r»i/2 point for Aa = —100 
GeV, tan /3 = 10, and sgn /i > 0. 




Figure 6: Relative error (<ti, /Lo) in the determination of the LSP decay length as a function of the mo CS> m\i% for Aq = —100 
GeV, tan/3 = 10, sgn /j, > 0, and an integrated luminosity of 100 fb _1 . The left (right) panel assumes no error (10% error) in 
the measurement of distance traveled by the LSP. 
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1 Y eqq' 




1 * Tqq' 


N eTV 




N TTV 


0.291 


0.106 


0.011 


0.087 


0.126 


0.061 



Table I: Reconstruction efficiencies for neutralino LSP decays for our benchmark point. For the r lepton only hadronic final 
states have been considered while the r decays into electrons and muons were included in the first two entries. 
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Figure 7: Expected error on the Br(x? — > £W + vZ) as a function mo ® mi/j for an integrated luminosity of 100 fb 



final state and simulated 100 fb -1 of data for several points in the mo <8> m.1/2 plane. As one can see, this branching 
ratio can be well determined in the regions of large production cross section, i.e. small m and m 1 / 2 - Although for 
heavier neutralinos the precision diminishes, still this branching ratio can be determined to within 20% in a large 
portion of the parameter space. In order to study the possibility of LHC to probe the atmospheric mass, we have 
evaluated ~Bt(x\ — > W£) + Br(x? — > Zv) appearing in Eq. [5] The W£ channel is obtained by first reconstructing 
displaced vertices with hadronic W decays, jj£, in the final state. Beside the cuts described in sections IlIII and ITVl we 
have applied an invariant mass cut on the jet pair: \Mw — \ < 20 GeV to disentangle the W-contribution to this 
final state. Afterward we get the branching ratio for W£ using 



Br(x? -> W£) 



Br(x? ->■ jj£) 



Ni 



</</' 



Br(W -> £v) 
Br(VF ^ qqt) 



(11) 



The Zv channel was calculated similarly by reconstructing the displaced vertices with hadronic Z decays, jjv, in the 



final state and properly rescaling it. Also here we have applied an invariant mass cut on the jet pair: 
GeV. 



\Mz-Mn \ < 20 



VII. LSP PROPERTIES AND ATMOSPHERIC NEUTRINO OSCILLATIONS 



As seen in Section II, the MSSM augmented with bilinear i?-parity violation exhibits correlations between LSP decay 
properties and the neutrino oscillation parameters 3l|, l32| , which are by now well measured in neutrino oscillation 
experiments 39j . In particular the squared mass difference Am§ 2 is connected to the ratio R32 between the LSP decay 
length and its branching ratio into £W and vZ; see the right panel of figure [5] In figure [5] we display the expected 
accuracy on the ratio R32 as a function of mo ® mi/j for an integrated luminosity of 100 fb -1 and assuming 10% 
precision in the determination of the LSP traveled distance. As we can see R32 can be determined with a precision 
20-30% in a large fraction of the mo^m^ plane and, as expected, the precision is lost for heavy LSPs. For small LSP 
masses the error on R32 is dominated by the uncertainty on the decay length, while for heavier LSPs the dominant 
contribution comes from the branching ratio determination due to the limited statistics. 
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Figure 8: Expected accuracy on the ratio R32 as a function of mo ® mi/2 for an integrated luminosity of 100 fb 

It is interesting to notice from the right panel of figure [5] that a measurement of R32 with 20-30% precision it is 
enough to determine the correct magnitude of Am| 2 using the BRPV-mSUGRA framework. Nevertheless, a much 
higher precision is needed to obtain uncertainties similar to the neutrino experiments such as MINOS/T2K [39! ]. On 
the other hand, the relation between the atmospheric mixing angle and the ratio of the LSP branching ratios into tW 
and fiW can lead to more stringent tests of the BRPV-mSUGRA model. In Ref. j^] it was shown that this ratio 
can be determined at the LHC with a precision better than 20% in a large fraction of the mo <8> mi/a plane. From 
figure [5] we can see that this precision is enough to have a determination for tan 2 #23 with an error similar to the low 
energy neutrino oscillation measurements. Looking from a different point of view, the collider data can be combined 
with neutrino data to determine the underlyin g pa rameters of the model. In this case collider and neutrino data give 
'orthogonal' information as has been shown in |46| . 



VIII. CONCLUSIONS 

We have analyzed the LHC potential to determine the LSP properties, such as mass, lifetime and branching ratios, 
within minimal supcrgravity with bilinear i?-parity violation. We saw that the LSP mass determination is rather 
precise, while the LSP lifetime and branching ratios can be determined with a 20% error in a large fraction of the 
parameter space. This is enough to allow for qualitative test of the BRPV-mSUGRA model using the i?32-Am 2 2 
correlation. On the other hand, semi-lcptonic LSP decays to muons and taus correlate extremely well with neutrino 
oscillation measurements of #23- 

In the BRPV model for low values of M1/2 one can have sizeable branching ratios into the final states erv and 
firis. These decays are potentially interesting for testing another aspects of the model associated with solar neutrino 
physics. As shown in [32| in regions of parameter space where the scalar taus are not very heavy, usually the loop 
with taus-staus in the diagram dominates the 1-loop neutrino mass. In this case the solar angle is predicted to be 
proportional to (ei/?2) 2 °c tan 2 6*©. Here, e = Vj ,tree e : with Vj' tree being the matrix which diagonalizes the tree-level 
neutrino mass. Note that vj ,tree is entirely determined in terms of the A^. In the BRPV model, RPV couplings of the 
scalar tau are proportional to the supcrpotcntial parameters Ratios of the decays Br(x" — > erv) /Br(x° ~ ^ f 1Tl/ ) 
are then given, to a very good approximation by Br(xi — > stv) fQi(x\ — > M ri/ ) x ( £ i/ £ 2) 2 - If the A4 where known, 
this could be turned into a test of the prediction for the solar angle. Note that in the limit where the reactor angle 
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is exactly zero and the atmospheric angle exactly maximal one obtains (ei/e^) 2 = 2(ei/e2) 2 . However, the Aj are 
currently not well fixed, due to the comparatively large uncertainty in the atmospheric angle. Thus the correlation 
between three-body leptonic decays of the ncutralino with tau final states and the solar angle has a rather large 
uncertainty. This prevents a stringent consistency test of the model using these decays. 

All in all we have shown that neutralino decays can be used to extract some of their properties rather well in 
models with bilinear i?-parity violation. Properties such as the decay length and the ratio of semi-leptonic decay 
branching ratios to muons and taus correlate rather well with atmospheric neutrino oscillation parameters. These 
features should also apply to schemes where the gravitino is the LSP and the neutralino is the next to lightest SUSY 



particle [47J, |48(. For gravitino masses in the allowed range where it plays the role of cold dark matter, its R-parity 
conserving decays are negligible compared to its R parity violating decays. The latter follow the same patter studied 
in the present paper, so that the results derived here should also hold. 
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